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Abstract 

Single electron tunneling into small superconducting islands is sensitive to 
the gap energy of the excitations created in the process and, hence, depends 
on the parity of the electron number in the island. We study these effects 
by analyzing the kinetics of the system. Since the interplay of Coulomb 
blockade and parity effects leads to a blocking of single electron tunneling, 
higher order tunneling processes become important. This is well-established 
for two-electron tunneling. Here we study processes of third and fourth order. 
We estimate the order of magnitude for these processes and discuss their 
relevance for recent experiments. 

PACS numbers: 73.40.Gk, 74.50.-hr,73.40.Rw 
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I. INTRODUCTION 



Experiments with systems including small superconducting islands showed effects which 
depend on the electron number parity in the island [0-0]. The effect has been interpreted 
as the even-odd asymmetry predicted by Averin and Nazarov |Q. It arises since single 
electron tunneling from the ground state, where all electrons near the Fermi surface of the 
superconducting island are paired, leads to a state with one extra electron - the "odd" one 
- in an excited state. In a small island, where charging effects prevent further tunneling, 
the odd electron does not find another excitation for recombination. Hence the energy of 
this state stays (at least metastable) above that of the equivalent normal system by the gap 
energy. Only at larger gate voltages another electron can enter the island, and the system 
can relax to the ground state. 

At higher temperatures, above a crossover value Tcr, only the e-periodic behavior typical 
for normal metal systems has been observed. This crossover has been explained by free 
energy arguments in Ref. [^. We had developed an alternative approach by studying the 
kinetics of the system Thus, we can describe nonequilibrium situations and derive 
the I — V characteristics of driven systems such as the SET transistors. We analyze the 
rate of tunneling of electrons between the lead and the island (reviewed in Section II) 
paying particular attention to the tunneling of the "odd" electron (or the tunneling of an 
electron into the partner state of the odd electron, resulting in an instantaneous condensation 
into a Cooper pair). The well studied single electron tunneling with rate F between a 
superconducting island and lead electrodes depends on the difference in the energy before 
and after the process. This includes the charging energy, but also the excitation energy Aj 
in the island. At low temperature processes which cost energy are suppressed. This renders 
the tunneling of the excited "odd" electron important. The tunneling rate 7 of this single 
electron is smaller by a factor l/N^^ff than the rate F. Here N^ff is the effective number of 
quasiparticle states available for an excitation at low temperature; in mesoscopic islands it 
is typically of the order of 10^. On the other hand, in an important range of parameters the 
rate 7 is not exponentially suppressed, since the excitation energy in the island is regained 
in this tunneling process. Hence 7 ~ Te^^^^'^'^ /N(,ff. Parity effects are observable as long 
as this single electron tunneling rate is relevant 7 ~ F, from which we obtain directly the 
cross-over temperature ksTcr ~ Aj/lnA^^e//. 

Our analysis of the tunneling rates can be applied to derive the I — V characteristics 
of normal-superconducting-normal (NSN) transistors (Section HI). In NSN transistors we 
combine single electron tunneling and 2e-tunneling via Andreev reflection to obtain a richly 
structured I — V characteristic, which reproduces quantitatively many features observed in 
the experiments of Refs. P,§J3l- I — V characteristic shows wide regions where both 
single and two-electron tunneling are blocked due to a combination of the superconducting 
gap and the Coulomb blockade, even at comparatively high transport voltages. The earlier 
experiments ||^ showed more structure, which could be identified as arising from the influence 
of the electromagnetic environment and could be removed by improved shielding [0. Even 
the latest experiments, however, show more structure than the theoretical results which 
account for single and two-electron tunneling. In particular the experiments show a ridge- 
like structure with a period e in the gate charge Qg at voltages eUtr ~ Aj. Cotunneling 
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processes are expected to produce an unstructured background rather than an effect which 
depends in a pronounced way on the gate voltage. Hence, the question arises whether higher 
than second order processes noticeably contribute to the characteristic. Here we consider 
coherent three-electron tunneling (Section IV). The process which turns out to be important 
is a coherent combination of a single electron transition and an Andreev reflection. The 
threshold voltage for this process and the order of magnitude of the rate match well with 
experimental results. 

Another interesting process is cotunneling of two electrons. As this process is of fourth 
order in the junction conductance its contribution is usually small compared to lower order 
processes. However, in the limit of small Utr both single electron tunneling and Andreev 
reflection (as well as their coherent combination) are suppressed due to the superconducting 
gap and Coulomb blockade. In this case cotunneling processes play the leading role. One of 
these processes - single electron cotunneling - has been investigated in Ref. P]. The process 
of two-electron cotunneling, i.e. the process of two coherent Andreev reflection events will 
be discussed in Section V. We will show that in the limit of low T and Utr either one- or two- 
electron cotunneling processes dominate depending on the system parameters. A discussion 
of our results is presented in Section VI. 



II. SINGLE ELECTRON TUNNELING RATES 

In order to review some basic results we consider an electron box assuming Aj < Eq- 
The results are easily extended to cover transport in SET transistors and systems with 
Aj > Ec- The electron box is a closed circuit consisting of a small superconducting island 
with total capacitance C = Cj + Cg connected by a tunnel junction with capacitance Cj to 
a lead electrode and by a capacitance Cg to a voltage source Ug. The charging energy of the 
system depends on Qg = Cgllg and the number n of charges on the island 

Ech[n,Qg)- — . (Ij 

The normal state conductance of the junction can be expressed by the tunneling matrix 
element and the normal density of states Ni/i{0) and volume Vi/i of the island and lead: 

1/Rt = 47Te^N,{0)VMm\T\Vn . 

The transition rate F"*" for a single electron tunneling (SET) process from the lead to the 
island, where the number of excess charges on the island changes from n to n + 1, is p,^] 

^^(Qs) = ^ / deM{e)fim ' ' ^ ' ^E^hiQ^)) • (2) 

6 Ihf J — CXD J — oo 

The 5-function expressing energy conservation depends on the change of the charging energy 
before and after the process 

SEchiQg) = Ech{n + 1, Qg) - Ech{n, Qg). 

If the distribution functions of lead and island are equilibrium Fermi functions the ex- 
pression for the rate reduces to 
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r+(Q,) = -h{SE,^{Q,)) \ . (3) 

The SET tunneling rate depends on the difference in the charging energy. In the super- 
conducting state it depends further on the energy gap, which enters via the BCS densities 
of states ^(e) = 6(|e| — Aj)|e|/i/e^ — A? of the island into the well known quasiparticle 
tunneling characteristic hieV) [p!0| . 

As long as the distributions are equilibrium Fermi functions the rate for the reverse 
process is given by the same expression as (0), however the sign of SE^h is reversed. Both 
satisfy the condition of detailed balance T~{Qg) = T~^{Qg)e^^'=''^''^'^ . At low temperature in 
the superconducting state the rates F^ are large only if the gain in charging energy exceeds 
the sum of the energies of the excitation created in the island ±6Ech + Aj < 0. They are 
exponentially suppressed otherwise. 

The assumption of equilibrium Fermi distributions is sufficient as long as we start from 
the even state. For definiteness let us assume that we started from n = and that the gate 
voltage is chosen such that \Qg\ < e. Hence, the rate of tunneling from an even to an odd 
state is 

F-(Q,) = F+(g,) . (4) 

However, in a superconductor with an "odd" unpaired electron, occupying a quasiparticle 
state above the gap, the distribution differs from an equilibrium one. Also the odd electron 
can tunnel back to the lead, enhancing the tunneling from odd to even states. Its initial 
energy in the island is at least Aj. This makes its tunneling rate large in a range of gate 
voltages where the competing processes, the tunneling of all the other electrons described by 
F~, are still exponentially suppressed. At finite temperature it is reasonable to assume that 
the odd state of the island is described by a thermal Fermi distribution but with a shifted 
chemical potential fsfii^i) = [e^^~^'^^^'^ + 1]^^. The shift in chemical potential is fixed by the 
constraint 

/oo 
deAf,{e)[fs^,{e) - fo{e)] . (5) 
-oo 

This reduces at low temperatures to 

6fi = \-T\nN,ff, (6) 

where 



N,ff{T) = N,{0)V,p7rA,T (7) 

is the number of states available for quasiparticles near the gap [0]. The tunneling rate from 
the odd state to the even state F"^ is given by the expression (0), however the island distri- 
bution function is replaced by /^^(e). For the following discussion it is useful to decompose 
this rate as 

F-(Q,)=F-(g,)+7(g,), (8) 
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where T has been defined above, and 

1 ^OO POO 

7(Q.) = ^ / d^lM{e^)[fs,{e^) - /o(e.)][l - foimi^i " - . (9) 

6 itf J —oo J—oo 

In the important range of parameters the rate 7 reduces to 

= 2e2/?,iV.(0)r. i for - - A. > T 

In comparison to F" the odd electron tunneling rate 7 contains a small prefactor l/Neff{T). 
On the other hand, there exists a range of gate voltages at low temperature where 7 is larger 
than r~, since it is not exponentially suppressed. The reason is that the gap energy of the 
island is regained in the process described by 7. 

For exp(— Aj/T) ^ 1 the ratio of the rates for two transitions is 

r^vr'^" = exp[(5E,^ + 5^l)/kBT] . (ii) 

The rates F*^" and obey a detailed balance relation, however, depending on the free 
energy difference, which in addition to charging energy involves the shift of the chemical 
potential 5/1 as well. This free energy difference coincides with that introduced in Ref. 

Above we described the range < < e where tunneling occurs between the island 
states n = and n = 1. The range e < Qg < 2e can be treated analogously. The 
tunneling now connects the states n = 1 and n = 2. In this case, except for the single 
electron tunneling F^, which creates further excitations, one electron can tunnel into one 
specific state {—k, — cr), the partner state of the excitation {k,a) which is already present. 
Both condense immediately; the state with two excitations only exists virtually. The latter 
process is described again by 'j{Qg)- The symmetry implies T'^°^°'^{Qg) = T^°/°^{2e — Qg). 
Since the properties of the system are 2e-periodic in Qg we have provided the description 
for all values of the gate voltage. 

Given the rates F^° and F°^ we can analyze where the transition between the even and the 
odd state occurs. We can describe the incoherent sequential tunneling of charges between 
the island and the lead by a master equation for the occupation probabilities of the even 
and odd states We{Qg) and Wo{Qg)- It is 

^^^^ = -F-(gjw^e(Q,) + F-(Q,)iy„(Q,) (12) 

with We{Qg) + Wo{Qg) = 1. The equilibrium solutions are 

w^e(o)(g,) = F-(-)(g,)/Fs(g,), (13) 

where Tj:{Qg) = T°^{Qg) + T^°{Qg). For F"'^ > T^" we have We{Qg) ^ 1, i.e. the system 
occupies the even state, while for F^° ^ F"*^ the island is in the odd state. Finally, the 
crossover temperature T^^ is determined by 7(e/2) ^ r[6E^{e/2)]). The result 

T,r = Ai/\nNeff{T,r) , (14) 

coincides with that found in Refs. Tcr can be interpreted as the temperature, at which 
the average number of excitations in the island is at least one. This result is valid also for 
SET transistors. 
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III. I -V CHARACTERISTICS OF SET TRANSISTORS 



The analysis presented above can be extended in a straightforward way to describe even- 
odd effects in SET transistors, which consist of an island coupled to two tunnel junctions (left 
and right) and a gate capacitance. In this system the charging energy depends on the gate 
voltage Ug and the transport voltage. The total capacitance of the island C = Cg + Ci + Cr 
defines the energy scale Ec = e^/2C. The energy differences for tunneling processes onto 
the island in the left and right junctions are 



SE^,,,^E,- <^'%^""\ (15) 



For clarity we assume a symmetric bias Ui = —Ur = Utr and Ci = Cr- We further define 

Qtr = CUtr- 

We focus our attention on NSN transistors with an energy gap which is larger than the 
charging energy Aj > Ec- In these systems the odd states have a large energy. Hence a 
mechanism which transfers 2 electrons between the normal metal and the superconductor 
becomes important. The Andreev reflection provides such a mechanism |Tl[]. The master 



equation description can be generalized to include also this process. In the limit considered 
the rate for Andreev reflection is given by the same expression as that for single charge 
tunneling (j^) with the following modifications [jl2| : (i) the charge transferred in an Andreev 
reflection is 2e, and the charging energy changes accordingly, (ii) the expression for hiV) 
is linear as for normal state tunneling, but (iii) the effective conductance is of second order 
= Rk I {4:NchRt)- Here Rk = h/e^ is the resistance quantum, and the number of channels 
Nch ~ 10'^ depends on the correlations of the electron propagation in the lead and the island 
13| , p!^ . An important conclusion is that Andreev reflection is also subject to Coulomb 



blockade |T2|]. Because of the similarity of the rate to that of single electron tunneling it is 



clear that the shape of the I — V characteristic due to Andreev reflection also takes a similar 
form as in normal transistors. At low temperatures a set of parabolic current peaks is found 
centered around the degeneracy points Qg = ±e, ±3e, . . . [|lT| 



I^{SQM = G\u,r-A^-^) . (16) 

Here 5Qg = Qg — e for Qg close to e and similar at the other degeneracy points. 

Depending on the gate and transport voltage we find different mechanisms to dominate 
the I — V characteristic. Upon increasing Utr we enter a regime where single electron 



tunneling dominates. This "poisons" or blocks the condition for Andreev scattering |]TT 
Hence after an initial increase with increasing Utr the current suddenly drops at the threshold 
voltage 

eQg 



eUtr = 2[Ec-^ + A, j (17) 

to the value limited by the odd electron tunneling, i. e. an "escape" current of the order 

1 

^^'^ ^ 2eRtN,{m ■ ^^^^ 
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For a detailed comparison of our calculations with available experimental data from Refs. 
we solved the master equation using the parameters of the experiments. 

Fig. 1 shows the I — V characteristic as a function of both gate and transport voltage. 
At small transport voltage we find 2e-periodic peaks due to Andreev reflection; the peaks 
at larger transport voltages arise due to subsequent incoherent steps of single electron tun- 
neling and Andreev reflection processes (which we will refer to as AQP cycles 0). As can 
be seen in Fig.l there are regions of Qg where even for transport voltages Utr ~ Aj/e both 
single electron and two-electron tunneling is blocked. In the experiment, however, a current 
is found in these regions. Especially, one finds ridge-like structures at Qg = 0, ±1, . . . which 
start at a transport voltage of Utr ~ 250fieV. We argue that these structures cannot be 
explained by the presence of several unpaired quasiparticles or inelastic cotunnehng. Con- 
sider the system at Qg = 0. There is no process of first or second order to alter the number 
of electrons in the island, hence, it cannot be explained why unpaired quasiparticles should 
appear on the island at the considered transport voltages. Moreover, increasing the gate 
voltage one approaches the AQP peaks and the probability to have unpaired quasiparticles 
increases, so one expects that a correspondingly large background has to be added to the 
current due to AQP processes. The experiment, however, indicates that instead of such a 
background a well-defined (e-periodic) structure is added to the features in Fig. 1. This 
motivated us to investigate higher order processes. 



IV. COHERENT THREE-ELECTRON TUNNELING 

An argument in selecting the important third-order contributions is that the processes 
may neither produce a large number of excitations nor change the number of electrons on the 
island considerably. This favors processes which are coherent combinations of two-electron 
tunneling in one junction and quasiparticle tunneling in the other (see the discussion below). 

The calculation of the rate is analogous to the calculation of the two-electron tunneling 
rate [jlT|. The rate for a process with two-electron tunneling in the left junction and one 



quasiparticle tunneling in the right can be expressed as 

li 



r(3) = |[ ^ \Mu.,r?f{il)f{iv){l - /(6fc))(l - fi^r) X 
ll'kr 

x6{^r + ek-^i~ 6' - ^eUtr + 6E,h) • (19) 



It contains the Fermi functions for electrons with energies C,i, C,i' in the left lead, C,r in the 
right lead and for excitations with energy on the island. We note the important fact that 
the process has a threshold voltage 

eUthr = '^{A, + 6E,h) . (20) 

Below this threshold the rate is exponentially suppressed. The matrix element Muik^ has 
the form 
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U'kr 



E 

k' 



rp(l) rp{l)rp{r) 

J-l-k'-'-lk J-rk V-k'Uk'Vk 



>fc + Cfc' + - 6 - eUtr)iek + ^r + 5^ch(-e) - eUtr/2) 

+ 5 other terms 



(21) 



where Uk,Vk denote the coherence factors in the island and T*^'/^) the tunnehng matrix 
elements of the left and right junctions. The combination of energy denominators written 
explicitly corresponds to quasiparticle tunneling in the first and Andreev tunneling in the 
second and third coherent step. We assume Utr ~ Uthr and therefore set ^i,^i',^r ~ and 
ek ~ Aj. Since we expect the rate to be important at integer values of Qg/e, we replace SEch 
by Ec- Carrying out the integrations and estimating the expressions for the parameters of 
the experiment |^ we find 



M, 



U'kr 



E, 



c 



(22) 



where the bracket {■ ■ ■)p denotes averaging over the directions of p. 

In the next step, we calculate the voltage dependence of the rate, assuming that T is 
much less than any other energy difference. The result is 



r(3) 



1 



leUtr - 5Ech 



A, 



elltr — SE, 



+ 2 



^eUtr-SE^h 



ch 



arcosh ( — 



A,- 



eUtr - 6E 



- 1 



ch 



A,- 



where 



An 



V = Yi\Mii'krnn'krNiiOyVi'NriO)VrN,{0)V,Ai 



(23) 



(24) 



The momentum average {Mii'kr)ii'kr is proportional to the product of the quasiparticle con- 
ductance of the right junction and the Andreev conductance of the left junction. Using the 
experimental parameters, without further assumptions concerning the averages of the tun- 
neling matrix elements, we estimate r] ~ 10^s~^. This compares well with the experiments. 

The master equation is completed with the rates F*-^-' for the various processes which 
change the island charge by ±e. The result is shown in Fig. 2. As we have expected there 
are no longer regions with blockade of the current at higher transport voltages. We also 
note that the ratio of the heights of Andreev and AQP peaks is changed. The simulation, 
however, does not reproduce the peaked structure of the experimental results, especially the 
ridges at = 0, ±1, . . . are missing. 

Finally, we return to the discussion of the criteria to select the important three-electron 
contribution. Other possible realizations are combinations of single electron tunneling and 
Andreev reflection in the same junction, coherent tunneling of three electrons with the 
creation of three quasiparticles in the island (tunnel events at the same or at different 
junctions) and combinations of elastic cotunneling and single electron tunneling. The first 
and the second of these realizations can be excluded, since they yield much higher threshold 
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voltages than given in (]20|). Furthermore, the experiment shows that elastic cotunneling 
plays no role for the I — V characteristic at transport voltages elltr ~ Aj. This is confirmed 
by an estimate using the theory of Ref. Hence, higher order processes including elastic 
cotunneling can be neglected. 



V. TWO-ELECTRON COTUNNELING 

As we have discussed above, the 2e-tunneling process through NS interfaces via Andreev 
reflection is to a large extent similar to the process of single electron tunneling between two 
normal metals. We can proceed further with this analogy and study the process of two- 
electron cotunneling in NSN transistors (similarly to the process of single electron cotunnel- 
ing in normal tunnel junctions |T^-|l8|). Let us consider two successive Andreev reflection 



events at two NS boundaries of the NSN transistor. One of these events corresponds to 
tunneling of two electrons from a normal metal into a superconducting island through one 
of the NS boundaries (with conversion of two normal electrons into a Cooper pair). For 
small values of the external voltage and at low T this process is energetically suppressed 
because it increases the charging energy of the superconducting island by ~ 4:Ec. The is- 
land, however, can be immediately discharged due to another coherent tunneling event of a 
Cooper pair from the island into the normal metal through the second NS boundary. Thus 
Coulomb interaction induces time correlation between two acts of Andreev reflection at two 
NS interfaces and the Coulomb blockade is lifted for such a process similarly to the case of 
single electron cotunneling. For a nonzero transport voltage Utr the above process provides 
a finite contribution to the current through the system. This contribution is proportional 
to but nevertheless it can become significant in the limit of small transport and gate 
voltages, i.e. in the case of Coulomb blockade of Andreev reflection at each of the two NS 
interfaces. 

In order to evaluate the two-electron cotunneling current through NSN transistors we 



shall make use of the method developed in Ref. [1/] for the case of single electron cotunneling 
in chains of normal tunnel junctions. It enables us to evaluate the contribution from the 
cotunneling process to the free energy Fg"* and then to calculate the current by means of 
the formula J|g* = 4eImF2™*. As the two electron cotunneling mechanism is important for 
small transport and gate voltages we restrict ourselves to the case Qg = and elltr <^ Eq. 
Proceeding perturbatively in 2 = ^q^i,2 i^i,2 the Andreev conductances of two NS 
junctions) and assuming that do not depend on the external voltage in the limit of small 
Utr one can easily show that Jg"* oc Ufj. in complete analogy with single electron cotunneling 
through normal metallic grains [P^5|-[T^. A rigorous calculation yields 



,eo. _ i^eraW^Ufr /arctan./5H^' 



(25) 



This result is valid provided the Andreev conductance of each NS interface is Ohmic. As 
it was demonstrated in Refs. [0,0 under certain physical conditions due to the proximity 
effect the Andreev conductance of NS interfaces may become non-Ohmic °^ ^tr^ the 
limit of small voltages. In this case we find /|g* = G^^^Utr, where [|T1| 
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i.e. the two-electron cotunneling conductance G^^^ turns out to be Ohmic. For reasonable 
experimental parameters we can estimate G^^^ ~ 10-^° Q-\ Thus the effect is in the 
measurable range and can be important in the limit of small external voltages. 

The effect of environment on two-electron cotunneling can be treated in complete analogy 
with that for the case of single electron cotunneling In the case of Ohmic external 



impedance Rx the above results for the cotunneling conductance should be multiplied by a 
factor ~ (2ef/f,.i?^C)4'='^-/^. 

Finally let us point out that the effect of two-electron cotunneling also occurs in SNS 
structures. The corresponding cotunneling conductance of such structures is identical to 
that for NSN systems with equivalent parameters. 



VI. DISCUSSION 

In conclusion, we have developed a theory of parity effects in small superconducting 
islands by analyzing the rates of various tunneling processes. The tunneling of the single 
odd excitation created in earlier tunneling processes is crucial at low temperatures. By 
comparing its rate 7 to the tunneling rate F^ of all the other electrons we explain the 
transitions and the crossover conditions in an electron box. Our approach can easily be 
generalized to derive the I — V characteristics of more complicated systems, such as SET 
transistors, which allows us to explain numerous details observed in recent experiments. The 
agreement with experiment can be improved by including coherent third-order processes. 

We considered a superconducting island which always remains coupled to the leads. In 
this respect our model differs from that considered in Refs. |[T9|-|2T|. These authors study a 



system with fixed electron number parity, either in the even or in the odd state, while single 
electron tunneling processes which allow for transition between the two states are forbidden. 
This restricts the excitation spectrum, leading to unusual thermodynamic properties. If the 
island remains coupled to the lead the excited states alternate between even and odd parity, 
resulting in rather different thermodynamic properties. 

Our analysis yields a richly structured I — V curves for NSN transistors and allows us 
to distinguish different charge transfer mechanisms important for different values of the 
external voltage. In particular, we found that in the low temperature limit the following 
processes are important for relatively large gate and/or transport voltages Qg and Utr'- 

a) Single electron tunneling (influenced by parity effects) p 

b) Two-electron tunneling (Andreev reflection) and 

c) Coherent three-electron tunneling. 

We have shown that the correspondence between the simulations and the experimental 
results is improved if the theoretical model takes into account three-electron processes. This 
model, however, does not reproduce the ridges in the I — V characteristic in [0. On the 
other hand, we mention that the results of Ref. |H] also show no signs of these ridges. 
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For small values of the external voltage and low T the processes a) — c) are suppressed 
due to the superconducting gap and Coulomb blockade. Under these conditions the system 
conductance is determined by one of the cotunneling processes: 

d) Single electron cotunneling ^ or 



Two-electron cotunneling |T^ 



As the corresponding contributions to the system conductance depend on different pa- 
rameters (in the case d), the conductance is inversely proportional to RfNi{0)Ai 0, whereas 
the conductance e) is defined by ([25|), (^6])) each of them can dominate depending on par- 
ticular experimental situation. 
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Figure Captions: 

Fig. 1: The current I{Qg, Utr) through a NSN transistor. Only single and two-electron 
tunneling has been taken into account. The parameters are chosen to coincide with those 
of Ref. §, Fig. 3a {Ec = 99/ieV", = 245/iel/). 

Fig. 2: The current I{Qg, Utr) through a NSN transistor including three-electron tunneling. 
The parameters are chosen to coincide with those of Ref. [0 , Fig. 3a {Ec = 99 fj,eV, Aj = 
245/ieV^); for the prefactor of the three-electron rate we set ?7 = 2 ■ 10^s~^. 
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